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ABSTRACT 


Particle-and-Force models for high-speed-computer calculations of 
fluid dynamics problems have several conceptual, advantages over most pre¬ 
sently used techniques; in addition, test calculations with several ver¬ 
sions have indicated considerable promise for wide application., For these 
reasons, the analysis presented here was undertaken to demonstrate the 
basis for expecting success in one class of such models (the completely 
statistical class) and in addition to show how to choose the interpar¬ 
ticle force function to correspond to a given material equation of state. 
The method of analysis employs a many-particle distribution function equa¬ 
tion whose reduction is accomplished through the assumption of two-body 
forces and of strong correlations with the mean. Primary results in¬ 
clude correspondence equations between conservative force function and 
equation of state, and between dissipative force function and viscous 
stress. Solutions are given for the first type of correspondence equa¬ 
tion, and comparisons are made with previous heuristic studies. 



Introduction 

In the pioneering work of Pasta and Ulam (1), it was shown that many 
features of the dynamics of fluids can be represented by a calculation of 
the trajectories of interacting particles, each representing a macroscopic 
fluid element. They studied in particular the problem of the instability 
at the interface between a heavy fluid falling into a lighter one, with 
the entire configuration resolved by only 256 points. The fluids were 
forced to move adiabatically, representing an unspecified equation of 
state; nevertheless, the results are surprisingly realistic and suggest 
that an extremely useful computing method could be derived as a generali¬ 
zation of their techniques. 

One such generalization that has received considerable application is 
the Particle-in-Cell (PIC) method (2) in which the particle dynamics are 
controlled by calculations that refer to an Eulerian mesh of cells through 
which the particles move. Another is the Nearest-Neighbor method of 
Kolsky ( 3 ), an elaborate and promising technique for which actual computa¬ 
tions are yet to be performed. Finally, there is the Particle-and-Force 
(PAF) method (4,5) for which preliminary tests have indicated potential 
usefulness. 

Of these modified approaches, the PAF 'method is most like the original 
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procedure of Pasta and Ulam. More conservative is Kolsky• s technique 

which, in effect, employs many-body forces among the particles, and for * 

which he is able to give a direct, mechanical derivation of the relation¬ 
ship between the representation and the true fluid dynamics. The PIC 
method is likewise conservative, resembling the well-known Eulerian meth¬ 
ods more closely than any other particle method. It is, however, inter¬ 
mediate between the Kolsky and PAF methods in that the accuracy of results 
depends to some extent upon the proper statistical averaging of some of 
its local variables (6); convergence and accuracy analysis is only partly 
mechanical and requires certain statistical assumptions. At the other 
extreme, the PAF method and its ancestral forms are in a certain sense 
completely statistical. This has the effect of rendering them capable of 
analysis, however, and the main purpose of this paper is to show how the 
methods of kinetic theory can be applied as a first step* 

Dynamics of the Particles 

Consider a set of particles labeled with indices i,j =1 N, and 
characterized by coordinates r. and velocities u • For this discussion, 
they all have the same mass, m. The force exerted by particle #i onto par¬ 
ticle #j is ?. and is assumed to be a function only of the variables de- 
i j 

scribing those two particles. It is further assumed that the appropriate dy¬ 
namics of the model is described simply by the usual law of mechanics, 

du* 

_ _A _ y* (1) 

m dt ^ ^ij' V 

where the * excludes the term i = j and may also restrict the sum to 
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certain neighbors of #j. The dynamics cannot follow completely the clas¬ 
sical laws of particle mechanics, however, since the elements of fluid 
that the particles represent must have another degree of freedom to allow 
dissipative effects. It is not necessary to specify at this stage the 
precise form of this extra variable, except that we shall require it to 
be a function of the entropy of the element, only; it is here designated 

by cp . We shall, however, be required to derive a precise form for the 
J 

equation of variation of cp., and in preparation, it is assumed that the 

J 

force function can be separated into two parts. 


“ r i/ ( vV * = Vm + hr 


( 2 ) 


where 


r iJ= r i’ r j|' 


{ -4 -4 x 

? - ( i: r i } 
iJ " r u ’ 


and cp. . is the average of cp. and cp. 
iJ i j 


The g_ force is introduced only to produce dissipation, so that the 

— 

cp. values will change only in circumstances in which g.. is nonzero. The 
0 j* J 

dissipative force must therefore be velocity-dependent, must vanish if 

there are no velocity gradients, and indeed for translational invariance 

must depend upon velocity as a function of u*. - u*. only. 

l 0 

To complete the dynamical description of the model, it is convenient 
to invoke the basic conservation laws from the viewpoint of fluid 


tt 
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dynamics. Conservation of mass is automatic; the particle masses are 

constant in time, and in addition the change of mass in any fixed volume 
exactly equals the amount flowing over the bounding surface. For con¬ 
servation of momentum it can be shown that it is sufficient that 

F = - F , from which it follows from Eq. (2) that cp = <p.. and 
ij ji J 1 


The energy conservation equation marks the point of departure from 

ordinary particle mechanics and furnishes the required equation for the 

changes of cp.. Let K. s (m/2)u..u* be the kinetic energy of particle 
J J J J 

while J. is its internal energy. The basis for the energy discus- 
J 

sion is that the rate of change of energy of a particle should be given 


by the rate that the other particles do work on it. This work rate is 
in turn given by the product of force by velocity. Overall conservation 
will require antisymmetry of this work flux, so that the velocity which 
carries it must be the average from the two particles. Thus we write 



It follows, as required, that the energy of any subset of particles 
changes only through work done on it by external particles. 

Now, by definition 


(3) 


dK. du 

—i = mu., 
dt j dt 



which, in combination with Eqs. (2) and (3) > 


can be put into the form 
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(4) 




* 

f 


dr . 
_JaL 
ij dt 


i 


+ 



* 



i 



This in turn may be compared with the identity 


dJ. 

dt 




*■ "ij 


at 


-he 


dj. . dr. . d<J, . dcp. . 

U iA + U M , 

.dr . dt dcp dt _r 
ij ij 


in which J . is the total internal energy of the ij pair, half of which 
is assigned to each particle„ 

This comparison shows that the internal energy of a particle changes 
according to two processes, one having to do with coordinate variations 
and the other with cp variations. The first change is of the potential 
energy relative to surrounding particles, while the second, referring to 
entropy changes, suggests that cp^ can now be set equal to the entropy it¬ 
self. Using the first law of thermodynamics, we are led to the result 
dcp. - c—i * 

'j 


dcp j 1 V*- (- -M 

”^j dt 1 " 2 L S ij* U i " V' 


(5) 


The temperature, T., is considered to be a function of cp. and of the mean 
j 0 

interparticle spacing, which dependence need not be further discussed 
here. We see as a consequent restriction that because cp. must increase 

J 

monotonically in time, the dissipative part of the force must be such 

that i? .(u* - uM is positive definite. 

0 ij i j 

To proceed, now, to describe an actual computing method, it would be 
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, to set up a 


necessary to prescribe in detail the forms of f^ and g ij 
te chni que for finite-difference time advancement consistent with sta¬ 
bility and accuracy requirements, and to determine the appropriate man¬ 
ner by which initial and boundary conditions could be incorporated. An 
example of these, together with discussions of results and interpreta¬ 
tions, is given in reference 4. There, however, choice of the force 


function, f , was made on a very loose basis, whereas here it will be 
ij 

seen that the nature of this function can be related to expressions for 


the continuum equation of state in a more satisfactory manner by the sta¬ 
tistical analysis. The choice of g^ in reference 4 was likewise not 
based on rigorous arguments, but again the statistical treatment shows 
how the force can be related to viscous forces (real or "artificial"). 


Statistical Analysis of Particle Representations 

It is known from the start that a particle calculation with only two- 
body forces cannot exactly represent the dynamics of a fluid. In reality, 
the force exerted by one element onto another must depend upon the rela¬ 
tionship of the first element to all of its neighbors; the more it is 
compressed by closeness of those neighbors, the greater will be the force 
it exerts. In a PAF-like method, however, there is no provision for this 
effect to be accurately represented at every instant, (in contrast, the 
Kolsky method with its many-body forces can properly represent the effect.) 
The only hope, then, for PAF-like methods is that the representation is 
good in some statistical sense. If this can be demonstrated, the numerous 
advantages of a PAF-like method can then be exploited. 




The procedure for determining whether or not the statistical averag¬ 
ing does give a proper representation can he outlined as follows# The 
starting point is a generalized Liouville equation appropriate to the 
study of the probable distribution of a set of dissipative particles. 

This is then reduced by the well-known integration method to the one- 
particle distribution-function equation which involves integrals of the 
two-particle function. (This complete reduction can be carried through 
o nly as a result of the two-body force assumption.) The resulting equa¬ 
tion is exact, but a study of its moments, from which the fluid-dynamic 
analogy is observed, requires some simplifying assumptions. In essence, 
these assumptions are, first, that there are strong interparticle corre¬ 
lations and, second, that there are strong correlations between individ¬ 
ual characteristics and those of the local mean. Surely if the particles 
are to represent macroscopic elements of fluid, then these assumptions 
must be good. Our results then show that the use of these assumptions 
allows the force-function conditions to be derived such that the par¬ 
ticle representation of the fluid dynamics is indeed statistic a l l y valid. 
Stated in other words, this means that if a choice of technique can be 
found such that each particle follows closely the local mean of motion 
(with no interpenetrations, strong oscillations, etc.), and if the force 
functions are properly chosen, the results should statistically repre¬ 
sent the fluid dynamics as desired. The remainder of this section shows 
how this analysis can be accomplished in detail. 

The basic entity for the discussion is the distribution function for 



the N particles, z^Cr^ *°° ^ ••• C P-| •*° ‘Tj ^^)> which is the 

joint probability density at time t for finding the particles in a con¬ 
figuration characterized by the specified variables. Normalization is 
such that integration over the allowed ranges of all the particle var¬ 
iables gives NJ, independent of time. (Normalization was chosen to be 
NJ rather than N or unity in order to conveniently count the permutations 
of the identical particles.) From z^ there can be derived subset dis¬ 
tribution functions 




•. o u ,cp 
n 1 



1 

(N - n)i 


r> 






and in particular we shall be concerned with z^ and 

In the generalized r, u*, cp space through which the particles move, one 
can imagine a closed surface defined to move in such a way that every 
point follows a particle trajectory allowable by the dynamics. Such a 
surface must then contain a fixed number of particles, so that the time 
derivative of the integral of z^ over all particle variables within the 
ranges contained by such a surface must vanish. Since this must be true 
for any arbitrary initial choice of surface, it follows that the inte¬ 
grand obtained from the time different at ion must vanish, leading to the 
generalized liouville equation 



j=1 






* 
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» 


where the subscripts of the divergence operators give the coordinates of 

variation, and the total time derivatives are constrained to refer to 
trajectories allowed hy the dynamics. These trajectory constraints are 
just those derived in the preceding section, Eqs. (1) and (5), plus the 
velocity definition, vl , = dr*. /dt» With these, the Liouville equation 
becomes 



To reduce this equation to a usable form, integrate over the complete 
ranges of the variables of all the particles except #1 and impose the con 
dition that there must be no flux of probability across the extremities 
of those ranges. Then, after some manipulation, there is obtained the re 
duced liouville equation for the one-particle distribution function. 



( 6 ) 


The first of the two strong-correlation assumptions, which asserts 
that two neighboring particles must not differ by much in their variables, 
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can be expressed in the form 


Z 2 ( ^ >* 2 *^ 

s z ] (? 1 ,3;,<(> 1 ,t)p(? 2 ,t)a(p 12 ,<|> 12 ,r 12 )8(^ - - ? 12 )6(<P 2 - <P, - f> 12 > 


In this expression, p(r*,t) is the local number density of particles, de¬ 


fined by 


p(r,t) =J&uJ&y z^u/p/t). 


( 8 ) 


while <j(p 12 ,cPi 2 , r 12 ) is a radial function, yet to be specified, which van¬ 
ishes outside of the range of interparticle forces; p 12 denotes the aver¬ 
age density between the particles, z*, on the other hand, is that part of 
the two-particle function which is important only outside the range of 
interparticle forces; indeed, we shall assume that its product with the 
force f unc tion vanishes everywhere. The 5-functions are the usual Dirac 
entities, while ? 12 and p lg are small quantities representing the ex¬ 
pected mean differences in velocity and entropy between the particles. 

The first approximations for them would be 




, —7 —7 \ 

(r - r )-V_ 

1 x J 


(f> - ?,)-v 


U 


9 


V 


(9) 

( 10 ) 


where u anH 9 are bhe local means of bhese quanbibiesj bhab is^ 

r 

dXLj&qu z 1 (r,u,qp,t). 



( 11 ) 
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pep a J du J dq> cp z^r,VL,<p,t), ( 12 ) 

so that they both are functions of ?and t. This first approximation 
will turn out to be adequate for ^ but not for c? 12 ; the reason for this 
is explained in the text following Eq. (23). Apparent lack of interchange 
symmetry in Eq. (7) will be resolved by restrictions on the form of z ^. 
Substitution of Eq. (7) into Eq. (6) results in the equation 


+ u-V z + V .(Pz) + SI- = 0 


_ o 
- °' 


- y r y —*) 

in which the subscripts have been dropped from z ^, u^ cp 1 , and r^ and 

P 3 1 f 21 + (14) 

x &(u 2 - - a!| 2 )&(cp 2 - <P 1 - P 12 )> 

8 * 5tF7 %r ( "2 ' “l )e(? 2 ) ° (p 12' ,,> 12' r l2 ) (1 5 ) 

X 5(u 2 - ?, - a 12 )6(9 2 - <P, - P 12 )- 

Note that z 2 does not appear. 

The second of the two strong-correlation assumptions states that 
the instantaneous velocity and entropy of any particle differ inappre¬ 
ciably from their local means; mathematically, 

z.,(r,u,cp,t) = p(r,t) 6 (u - u)B(<p - 9 ). O 6 ) 

It can be shown that this choice for z 1 makes z 2 symmetric as required. 

It is useful to consider the normalization properties of z 2 and z^ 


-15- 



in the forms given by the strong correlation assumption. Thus, we must 
have 

JJJ z 2 dr 2 du 2 dcPg s (N - IJz^ 

.[[[ z i ^ d? i atp i = N< 

From Eq. (1 6 ) it is seen that the second normalization requirement 
becomes 

J p(r) dr = N, 

which only demonstrates consistency and gives no added information. With 
Eq. (7), the first normalization condition becomes 

1, Iff Z i ^2 *4 d<P 2 + / p(? 2> ° (p 12’”l + \ »12> r 12> ^2 * " - ’• 

Now z^ is considered to vanish within some region about particle #1, while 
-a vanishes outside that region. Furthermore, in the region in which z^ 
does not vanish, we may suppose that the two-particle correlation is weak, 
so that Zg = z^ (1) z 1 (2), which nomenclature signifies the product of the 
separate one-particle functions. Thus the normalization on z 2 becomes 

f p{r 2 ) d? 2 +Jp(? 2 ) »(P 12 ^i + I P 12' r 12 )d5 2 = N " 1# 

The primed integration is over the exterior region, and thus equals 
N - (N + 1), where N is the number of neighbors of particle #1 implied 
by the * sum of Eq. (1). It therefore follows that 



♦ 


4 
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To lowest order, then 

r 2 * 

if.it / r ps(p,cp,r) dr = N , 

J 0 

3 

which implies that a is a function of r p only, together with the condi¬ 
tion 

/ °(S) = N*. 

3 J o 

Here v is the cutoff value of pr heyond which a vanishes. Since the cor¬ 
responding radius must include N + 1 particles, we obtain for v the value 

v - & (N * + ')• 

This completes the preparation for the transition from particle dy¬ 
namics to fluid dynamics. The treatment resembles that used in the kinetic 
theory of gases, but differs in several essential respects. In particular, 
the two strong-correlation assumptions, both of which are appropriate to 
the present situation (and indeed are essential properties of the particle 
dynamics if the fluid-dynamics analogy is to be valid), will be seen to re¬ 
move from the statistics those fluctuations which in gas dynamics produce 
internal energy, heat transport, viscosity, diffusion, etc. These have, 
instead, been replaced by force functions which return the missing gas 
properties, but in a manner appropriate to the correlated motions of macro¬ 
scopic fluid elements. Still missing from the present treatment are the 
phenomena of heat conduction and diffusion. At least the former (and 
perhaps the latter) could be included in the present model by addition of 



appropriate interparticle fluxes* 

The next step, now, is to make the actual transition to fluid dynam¬ 
ics. This is accomplished hy taking appropriate moments of Eq. (13), in 
which z 1 is of the form given in Eq. (16). The simplest moment is ob¬ 
tained by integration of Eq. (13) over u* and cp. The result is the fluid- 
dynamics mass conservation equation: 


( p\x) = 0 . (17) 

r 

This result, which could as well have been obtained from Eq. ( 6 ), does not 
depend upon the strong-correlation assumptions. 

Next is the velocity moment which, after some manipulation, reduces 


to 


+ V ‘(pu3 = pP, 
ot 

r 


( 18 ) 


in which 1? means i? evaluated at u = n and 9 = cp. Here the consequences of 
the strong-correlation assumptions are clearly manifested in the absence 
of the gas-dynamics fluctuation pressure. 

Finally, the cp moment of Eq. (13) can be put into the form 


+ V_»(puqp) = pQ, (19) 

r 

where Q means Q evaluated at u = u and cp = cp. As expected, this shows 

that the entropy will change along the motion of a fluid element if and 

only if the dissipative force, g^j, is nonzero. 

■=♦ — 

The next step is the reduction of pP and pQ to more meaningful forms. 


9 
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from which, for example, the relationship between equation of state and 

force function can be derived. To proceed, it is necessary to examine 

the force functions in somewhat more detail. The nondissipative force 

has already been specified to be of the form f. . - f(r ,cp. ,)j in addi~ 

fj fj ij 

tion the dissipative force is now specified to be of the form 

g*. . = gj^r t -jT. .,u. - u*.), which is the most general form consistent 
lj i> J 

with symmetry and invariance requirements. 

Substitution into Eqs. (14) and (15) and insertion of mean values 
for u and cp result in the equations 


P m J dX 2 S 21 f ^ r l2 ,Cp + 2 ^12* + ^ r i2 ,S 21 ,CP + 2 P 12 ,a i2^ 


X p(r 2 ) o(p 12 ,cp + I P 12 »r 12 ), 

« ■ + i P 12'“l2> - “12 P(? 2 ) <l(p 12' <P + i VV’ (2,) 

In the integrand, those quantities not depending on the position of 
particle #2 are no longer subscripted, consistent with the dropping of 
subscripts in the moment equations. 

Consider, further, the integral in Eq. (20). It can be reduced 
through the expansions 


f( V* + 5 P 12> 


S(^ 2 ,t 21 , v + l^ 2 ,a 12 ) . 


p(? 2 ) 


x. . 1 df 

* + 5 e i2^ + “l2 - ( W -’S 

ocp a 

P + r 12 i 12-( V j) > 
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- 1 \ •\ * f— \ feo\ . 1 R do 

= C, + 2 r i2 S 12\ V -fAV + 2 P '2s’ 


in which the functions f , g^ and o on the right are evaluated at 

2 = ° f 2 = °* ^ *2 ~ ^ except that they remain of r i 2 

in their explicit dependence on that variable). The subsequent lengthy 

manip ulat ion consists of evaluating the angular parts of the integrations 

and the collecting of terms, all leading to the result 


p p = - V_ 


~ - & ■ f r 5 a[p(r),cp(r),r] f[r,cp(r)] dr 


M JO 


4^- “l2°( V -#/ dr 2* 


The second integral can be further reduced through the assumption that 

g = a 12 G(r 12 ,9) + s* 21 (s > 21 *« 12 ' ) H ^ r 12^ e ^ 

The first term is directed along the line of velocity difference while the 
second is along the interparticle line. The first term does not necessar¬ 
ily conserve angu la r momentum microscopically while the second one does. 
This expression for g is the most general reasonable combination to first 
order in velocity differences. We now let 


a = lr 12 s 12 -V 


\ -> 


— (r i* *V ) t 

2 V 12 12 V 


which contains the next step in the expansion beyond that in Eq. (9) , re¬ 
tained because all integrals vanish identically when only the first step 
is included. A further lengthy but straightforward set of angular inte- 
grations reduces the expression for pP to its final form 
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pP = - f x 5 0[x 5 p(r^] f[x,cp(r)] dx 


0 


(24) 


+ (A + p) vJ^V^u^ + |i vj 2 ^ 


r r 


in which 


d = 2, r| ' ^ > 7" x^a[x 5 p(i^)] (50 + H) dx, 
- iprtp 2 (?) f% r 3 i TI j 

>> + d = —T^iT J x a ^ x H dXe 


(25) 

(26) 


The normalization result has been incorporated into the form of cr, and x 
is a dummy integration variable. The brace factor in Eq. (24) is the 
pressure given by the equation of state as p(p,cp) > while p is the first 
viscosity coefficient and A + 2p/3 is the second. 

‘ft 

We have seen that o(|) =0 for | > v = (3/4n)(N +1); presumably 
o(g) -»0 as g -»0. We thus take, as a first approximation. 


o( |) = a = const for £ < g < v f 


= 0 otherwise, 
where o and £ are determined by 

7 fo ° d5 = '* 

T f* 035 - “*• 

This, then, gives a = 1, £ = 3/4jc. 
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We thus have the integral equation to solve: 


f x 5 f(x,cp) dx = p(p,q>), 

3m J y. 

X 1 

where x, = (3A*p ) l/5 and x 2 = [3(N* + 1>A«p] l/5 ; the result is the cor¬ 
respondence between force function and equation of state. Note that this 
equation applies to calculations in three space dimensions. An analogous 
derivation for two space dimensions results in a similar correspondence 
equation 

2 p X 2 

J x 2 f(x,cp) dx = p(p,cp), 

X 1 

1 /2 1 f 2 

in which we now have x 1 = (1 /pot) and x 2 = [(N + 1)/p«] 

Before discussing the solution of these integral equations, it is 

useful to consider two aspects of their derivation. The first concerns 

interpretation of N*, the number of neighbors of a particle. If these 

neighbors were chosen to include all within a certain fixed prescribed 

r adi us, then our assumption of two-body forces in the Liouville equation 

* 

reduction would indeed be valid. In that case, however, N would be a 
function of density, and the subsequent analysis would not be strictly 
correct. On the other hand, if N* were to include a fixed number of 
nearest neighbors, then the two-body force assumption would no longer be 
strictly valid, and the reduction itself would give only an approximation. 
In either interpretation, however, the inaccuracies must be small since 
the results as applied to polytropic gases have proved to be good, both 


'm 
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in conceptual form and in accuracy of application. 

The second aspect of the derivation concerns the meaning of cp, to 
which a correspondence with entropy has been assumed. A review of the de¬ 
rivations shows that this correspondence is required only in the analysis 
of the dissipative part of the interparticle forces. Insofar as the 
equation-of-state correspondence is concerned, as given in Eq. (27) or 
Eq. (28), cp could be any appropriate second variable in the equation of 
state. In particular, we shall find it convenient (and in some respects 
even necessary ) to interpret cp as the specific internal energy, the 
reason being that only with that association can energy be rigorously 
conserved in the finite-time-interval solution of the particle-dynamics 
equations. 

Consider, now, the solution of Eq. ( 27 ) for a polytropic gas. 

First, with cp = entropy, we have 

y C P/ C V 

P = Ap' e , 

in which c v is the specific heat and A is a constant. The solution for 
f is then 

2-d-y ^v 

f(x,cp) = Kx ' e > 


in which 


K = 


9mAf2 - y) 



3(N*+1) 


, 2-7 



This is a generalization of a result given previously (reference 7 > 
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page 37), which had been derived from single energy considerations. If, 
however, <p = specific internal energy, then 

P = (7 “ 1 ) P9, 
and the solution for f is 

f = 2E ( r - 1 ) $ , ( 29 ) 

N X 

in which q is the number of dimensions. [This is, in other words, the so¬ 
lution of both Eq. ( 27 ) for which q = 3 and Eq. (28) for which q = 2.] 
Comparison with other previous results (reference 4, page 25 ) shows exact 
agreement if one chooses N = 2q (as was indeed specified in that refer¬ 
ence, page 23 , for the primary number of neighbors). 

To find the general solution for Eq. (27) or Eq. (28), consider the 
problem of solving for $(x,cp) the equation 
pkr } 

/ ^(y^) ay = ^(0,9)• 

J a 

The two equations of interest are both special cases of this, and 
can be identified with it by the proper choices of <&, >|r, fl, and k. Note 
in particular that for Eq. (27), k = (K + 1) ' , while for Eq. ( 28 ), 
k = (N* + I) 1 / 2 . In both cases, therefore, k > 1. By differentiation, 
the equation becomes 

kO>(kfl) - <X>(ft) = ~ , 
from which it follows that 
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®(n) 




where n is any positive integer. In many circumstances, the equation of 
state will be such that we may let ty -♦ 08 and arrive at the solution 


<t(a,cp) = 


an 


00 

^*(nk' J ,q>). 


(30) 


j =1 


As an example, consider Eq. (28) with a polytropic gas in which cp is spe¬ 
cific internal energy. This leads to the identities 

n = (i/pn) 1 / 2 , 

\}r(n) = 2m(7 - 1 ) n 2 , 

®(y,<p) = y 2 f(y*<p)> 

k = (n* + i) 1//2 . 


and. the solution 

««*> - 4j • 

j=1 K 

*K* 

The sum converges to the value 1 /N so that the result agrees exactly with 
that of Eq. (29). 

In general, for two dimensions 

a - 0 /p*) 1/2 , 

\Kn,q>) = 2mjt n^p(—^,q>) > 

^ itn 

*(y,<p) = y 2 f(y/p ), 
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while for three dimensions. 


0 = (3/W) 1 ^, 

= | nmn 6 p (/ ^ 3^) > 

®(y,<p) = y 5 f(y>9), 

k = (N* + l)^ 5 . 

Note that in "both cases, a requirement for convergence of the sum in Eq. 
(50) is that p/p 2 -»0 as p -* », a requirement satisfied by all gases 
under ordinary circumstances. Materials whose equations of state do not 
satisfy this criterion will require alternative methods for solution of 
the correspondence equation. 

Conclusions 

It has been shown that if a Particle-and-Force computing method 
with two-body forces can accomplish strong correlations between particles 
pnri of each particle with the local mean, then the use of a force func¬ 
tion given by the correspondence equation can result in statistical repre' 
sentation of the true dynamics of a fluid. In addition, if some pre¬ 
scribed dissipative force is required to accomplish the correlations, 
then the analysis shows the relation of this to viscous stress (true or 
"artificial"); or conversely, if a certain viscous effect is desired, the 



corresponding dissipative force can be derived. 

Derivation of an actual working method requires considerably more 
than the statistical analysis presented here; some preliminary efforts 
have been described in the references, and at present Group T-3 of this 
Laboratory is engaged in a program of development which will be reported 
later. 
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